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We consider a two-dimensional fully frustrated Josephson-junction array, which is driven uniformly 
by oscillating currents. As the temperature is lowered, there emerges a dynamic phase transition 
to an ordered state with nonzero dynamic order parameter for small currents. The transition 
temperature decreases monotonically with the driving amplitude, approaching zero at a certain 
critical value of the amplitude. Above the critical value, the disordered phase and the dynamically 
ordered phase are observed to appear alternatively. The characteristic stochastic resonance behavior 
of the system is also examined, which reveals that the resonance behavior of odd and even harmonics 
can be different according to the zero-temperature state. 

PACS numbers: 74.50.+r, 74.25.Nf, 05.40.-a 
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I. INTRODUCTION 



Dynamic responses of cooperatively interacting many- 
body systems to time driven perturbations are ex- 
tremely important technologically and involve intrigu- 
ing physics, resulting in intensive investigation in recent 
years. Among the well-known systems is the simple-. 
netic Ising model under oscillating magnetic fields n 
It has been revealed that spontaneous symmetry break- 
ing Aftkes place at a finite strength of the oscillating 
fieldua and that in two dimensions the transition be- 
longs to the same universality class as the equilibrium 
two-dimensional (2D) Ising modelE The stochastic reso- 
nance (SR) phenomena have also been studied.^ 

Another interesting example is the dynamics of 
the fully frustrated Josephson-junction array (FFJJA), 
which possesses the same Z2 symmetry as the Ising model 
in addition to the continuous U(l) symmetry. Unlike the 
Ising model, the FFJJA has real intrinsic dynamics de- 
rived from the Josephson relations and thus grants direct 
experimental realizations. There have bosixsome studies 
on the dynamic properties of the FFJ JAoEl and recently 
dynamic transitions have also been investigated in the 
FFJJA driven by uniform dc currents^ or staggered os- 
cillating magnetic fields. El Similarly to the kinetic Ising 
model, the FFJJA driven by staggered oscillating fields 
has been shown to exhibit a dynamic phase transition and 
an SR behavior. In the presence of weak fields, the dy- 
namic phase transition belongs to the same universality 
class as the equilibrium Z2 transition in the fully frus- 
trated XY (FFXY) model. At strong fields, in contrast, 
a different universality class has been suggested. 

In this paper, we investigate the FFJJA in the pres- 
ence of uniform oscillating currents, which is easier to 
realize in experiment than the system under staggered 
oscillating magnetic fields. The dynamic properties of 
the system is examined with attention paid to the dy- 
namic transition and the SR phenomena. At low tem- 



peratures the chirality in the FFJJA displays antiferro- 
magnetic ordering. To describe the dynamic transition 
associated with the antiferromagnetic ordering of the chi- 
rality, we introduce the dynamic order parameter, which 
is the staggered magnetization averaged over one period 
of the uniform oscillating current, and investigate its be- 
havior with respect to the amplitude and the frequency 
of the driving current as well as to the temperature. At 
zero temperature, as the driving amplitude is raised, the 
dynamically ordered phase and the disordered one are 
observed to appear alternatively. As the temperature is 
increased, the ordered phase undergoes a dynamic phase 
transition to the disordered phase. Through the use of 
the finite-size scaling analysis, the dynamic transition is 
shown to belong to the same universality class as the 
equilibrium Z2 transition present in the FFXY. We also 
examine the SR phenomena in the system. As manifested 
by the signal-to-noise ratio (SNR) at several harmonics, 
odd and even harmonics are revealed to exhibit different 
SR behaviors. It is discussed in view of the variation 
of the zero-temperature states with the driving current 
amplitude. 

This paper consists of five sections: Section II in- 
troduces the fully frustrated Josephson-junction ar- 
ray, driven uniformly by alternating currents, together 
with the dynamic order parameter, describing the dy- 
namic transition of the system, and presents the zero- 
temperature behavior. In Sec. Ill the behavior of the 
dynamic order parameter at finite temperatures is inves- 
tigated in detail, on the basis of which the dynamic phase 
diagram is constructed on the plane of the temperature 
and the driving amplitude. The nature of the transi- 
tion is also examined and the equilibrium Z2 universality 
class is identified. Section IV is devoted to the power 
spectrum of the staggered magnetization and the corre- 
sponding signal-to-noise ratio, which reveals the charac- 
teristic resonance behavior of the system with the appro- 
priate values of the parameters. Finally, a brief summary 
is given in Sec. V. 
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II. FULLY FRUSTRATED ARRAY DRIVEN BY 
ALTERNATING CURRENT 

To begin with, wc consider the equations of motion for 
phase angles {4>i} of the superconducting order parame- 
ters in grains forming anLxL square lattice with unit 
lattice constant. Within the resistively-shunted-junctiop. 
model under the fluctuating twist boundary conditions,tl 
they read: 



dt 



sm 



A) + 



= 0, 



(1) 



where the primed summation runs over the nearest neigh- 
bors of grain i and the thermal noise current r/ij satisfies 
(%(*)'JhW) = 2TS(t—t')(8ik5ji—5aSjk) at temperature 
T. We have used the abbreviation 4>ij = 4>i ~ 4>j ~ 
and Tij = Ti — Fj with r, = (x i: yi) denoting the position 
of grain i. Note that for nearest neighboring grains 
is a unit vector since the lattice constant has been set 
equal to unity. We have also written the energy and the 
time in units of hi c /2e and h/2eRi c , respectively, with 
the critical current i c and the shunt resistance R. The 
dynamics of the twist variables A = (A x , A. y ), which are 
included in the fluctuating twist boundary conditions to 
allow fluctuations in the phase difference across the whole 
system, is governed by the equations 



dA x = _L_ 

dt ~ L 2 



sm((f>ij - Ajb) + rj A;i; - I sin(fii), 



dA y I 
~dT 



— ^ Sm(0y - Ay) + T] Ay , 

(#)«, 



(2) 



where Yluj) denotes the summation over all nearest 
neighboring pairs in the a (= x, y) direction, rjA a satisfies 
(VA a (t + ^)VA a ,(t)) = {2T/L 2 )5(T)5 a , a <, and the oscil- 
lating current Iq sin(rtt) is injected in the x direction. In 
the Landau gauge, the bond angle Aij, given by the line 
integral of the vector potential, takes the values 



A; 



for Tj = Ti + X 







nxi for r ; = Ti 



For the study of the transition associated with the Z2 
symmetry in the FFJJA, it is convenient to consider the 
chirality 



C(R,i) = sgn 



^sin^y(i) -Tij ■ A(t) 



and the staggered magnetization 



mi 



(3) 



(4) 



R 



where denotes the directional plaquette summation 
of links around the dual lattice site R = rj + (l/2)(x+y). 



To probe the dynamic transition in the presence of an 
oscillating uniform current, we define the dynamic order 
parameter as the staggered magnetization averaged over 
one period of the oscillating current 



Q = — 



m(t)dt 



(5) 



In the numerical calculation, the sets of the equations 
of motion in Eqs. (|l|) and (||) are integrated via the mod- 
ified Euler method with time steps of size At = 0.05. 
We have varied the step size, only to find no essential 
difference. Typically, data have been averaged over 5000 
driving periods after the data obtained from the initial 
1000 periods discarded; the appropriate stationarity has 
been verified. In addition, five independent runs have 
been performed, over which the average has also been 
taken, and systems of size up to L = 32 have been con- 
sidered. 

In the FFJJA, it is well known that there are two kinds 
of antiferromagnetic chirality ordering at zero tempera- 
ture. The injection of uniform oscillating currents plays 
the role of tilting sinusoidally the 2D "egg-carton" lattice 
pinning potential. Accordingly, if the driving amplitude 
Iq of the injected currents is large enough to overcome 
the lattice pinning potential, oscillations between the two 
ground states are expected to take place. In Fig. [I], we 
display the time evolution of the staggered magnetization 
m(t) at zero temperature, evolved from the initial con- 
dition m(i=0) = 1 for driving frequency Q/2tt = 0.08. 
For small Iq, the system is shown to stay in the state 
with m(t) = 1 in spite of the driving current. When the 
amplitude exceeds the value rj 1.01, the system is 
driven out of the m = 1 state and the staggered mag- 
netization oscillates between m = ±1; this arises since 
the chirality lattice moves collectively over the lattice 
potential under large current driving. In case that the 
driving induces the chirality lattice to move over only 
a single lattice barrier in the first half of the period, the 
system oscillates between the two ground states symmet- 
rically (see the staggered magnetization for Jo = 1.50 in 
Fig. ||). On the other hand, further increase of Iq al- 
lows the chirality lattice to go over another barrier and 
yields asymmetric oscillations of the system between the 
ground states, as demonstrated by the staggered magne- 
tization for Iq — 2.0 in Fig. |l|. The symmetry is restored 
and broken alternatively with the increase of Iq , which is 
manifested by the zero-temperature behavior of Q with 
Iq in Fig. ||. As expected in Fig. [j], the dynamic order 
parameter Q, which is unity below Iq^°\ is shown to de- 
crease abruptly to zero for I > I (0 \ Further increase 
of Iq reveals that the dynamically ordered state (Q > 0) 
and the disordered one (Q = 0) appear alternatively in 
the system. It is observed that the saturation value of Q 
is less than unity and that the peak value of Q in each or- 
dered region decreases monotonically with Iq. Note also 
that the increase of the driving frequency f2 enhances the 
value of Iq^ and moves the dynamically ordered regions 
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FIG. 1: Time evolution of the staggered magnetization at 
zero temperature for driving frequency Q./2n = 0.08 and am- 
plitudes Jo = 0.98, 1.03, 1.50, and 2.0 from above. For clarity, 
the data for lo = 0.98 and 1.03 have been shifted upward by 
3.5 and 2.5, respectively, whereas those for Iq — 2.0 shifted 
downward by 2.5. 



to larger values of Iq, which reflects that a higher fre- 
quency implies a shorter period and accordingly, larger 
driving is necessary to move the chirality lattice over the 
lattice barrier within the (shorter) period. Except for the 
scale, however, the overall behavior at zero temperature 
does not change qualitatively with fl. 



III. DYNAMIC TRANSITION 

We next investigate the effects of the temperature. At 
sufficiently high temperatures, thermal fluctuations are 
dominant so that the influence of the driving current and 
the lattice pinning potential can be neglected. Conse- 
quently, at high temperature, the staggered magnetiza- 
tion should fluctuate randomly with time, resulting in 
the vanishing dynamic order parameter. On the other 
hand, at low temperatures, where thermal fluctuations 
are small, the system is expected to be driven mainly 
by the oscillating uniform currents, exhibiting behavior 
similar to that at zero temperature. Accordingly, we ex- 
pect a phase transition between the dynamically ordered 
phase and the disordered one as the temperature is var- 
ied. Figure ^, in which the ensemble average of the dy- 
namic order parameter (Q) is plotted as a function of 
the temperature T for various driving amplitudes, ex- 
plicitly shows the existence of such a transition between 
the dynamically ordered phase and the disordered one. 
The dynamic order parameter, starting from zero at high 
temperatures, develops as the temperature is lowered. It 
then grows rapidly and saturates eventually to the zero- 
temperature value as the temperature approaches zero. 
In particular, for Iq above Jo j Fig- P)(b) shows that the 
zero-temperature value of {Q} is reduced rapidly with the 
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FIG. 2: Dynamic order parameter as a function of the driving 
amplitude at zero temperature in the system of size L = 16 
with (a) n/2-K = 0.08; (b) n/2iv = 0.16. 



increase of Iq. Here, it is interesting that the apparent 
peak of {Q) is observed at a certain finite temperature 
rather than at zero temperature. 

To estimate the ptransition temperature, we consider 
Binder's cumulantsu 



U L = 1 



(Q 4 



3<Q 2 } 2 ' 



(G) 



Since this quantity is size-independent at the transition 
temperature, the transition temperature can be deter- 
mined from the crossing point of Ul for several sizes L. 
In Fig. ||, we plot the resulting phase diagram on the 
T— Iq plane, displaying dynamic phase boundaries. The 
transition temperature initially decreases monotonically 
to zero as the driving amplitude Iq is increased. How- 
ever, further increase of Iq drives the system into an- 
other ordered region, where the transition temperature 
first grows with the driving current, then reduces to zero. 
The zero-temperature results shown in Fig. ^suggest that 
these additional ordered regions emerge repeatedly with 
the amplitude Iq raised although the peak values of the 
transition temperature in these regions should decrease 
with Iq. As the driving frequency fl is increased, on 
the other hand, the transition temperature also increases 
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FIG. 3: Dynamic order parameter as a function of the 
temperature in the system of size L = 16 for various driv- 
ing amplitudes (a) I = 0.3(D), 0.5(O), 0.8(A); (b) I = 
1.02(D), 1.04(O), 1.06(A). Error bars have been estimated 
from standard deviations; in (a) they are not larger than the 
sizes of the symbols. Lines are merely guides to eyes. 



and the dynamically ordered region expands on the T—Iq 
plane. It is thus concluded that high driving frequency 
helps to maintain the dynamic order, which is consistent 
with the zero-temperature result in Fig. ||. It should also 
be noted here that cooling and heating curves for the 
dynamic order parameter do not exhibit any appreciable 
hysteresis even in the strong-current regime. 

To probe the nature of the transition, we consider the 
scaling relation for the dynamic order parameter 



(Q) = L-V"f ((T - TJL 1 /") 



(7) 



and plot {Q)L / V versus (T — T c )L l / v for various sizes. 
Figure || exhibits the resulting scaling plots for the driv- 
ing amplitude (a) Iq = 0.3 and (b) Iq — 2.2, with the 
critical exponents: (a) v = 0.82 and f3/v = 0.11; (b) 
v = 0.95 and f3/v = 0.13. Since the critical exponents in 



FIG. 4: Dynamic phase diagram on the T—Iq plane for driv- 
ing frequencies Q./2n — 0.08(D) and 0.16(A). The boundaries 
are determined by the crossing points of Binder's cumulants 
for size L — 8, 16, and 24. Again error bars have been esti- 
mated from standard deviations and lines are merely guides 
to eyes. 



(a) agree well with those,£ojr the equilibrium Z 2 transi- 
tion in the FFXY model,EHL3 the nice collapse of the data 
suggests strongly that the dynamic transition in the first 
ordered region [i.e., Iq < Iq ^ in Fig. ^(a)] belongs to the 
same universality class as the equilibrium Z2 transition of 
the FFXY model. Similar conclusion was also reached in 
the FFJJA under weak staggered oscillating fields xi On 
the other hand, the values in (b), corresponding to the 
second ordered region, appear somewhat larger than the 
critical exponents for the equilibrium Z2 transition. How- 
ever, due to large fluctuations, the numerical accuracy in 

(b) is not so good, yielding reasonable collapsing behav- 
ior down to v = 0.87 and (3/v = 0.12. Further, note that 
there does exist rather large discrepance, among the equi- 
librium values reported in literature,EE30 which may be 
summed up as follows: (3 = 0.11 ±0.03 and v = 0.9 ±0.1. 
In view of this uncertainty, it is also possible to con- 
sider the values in (b) to be consistent with those for the 
same (equilibrium Z2) transition. Accordingly, we pre- 
sume that for all driving amplitudes the dynamic transi- 
tion in the system belongs to the same universality class 
as the equilibrium Z2 transition present in the FFXY" 
model. 



IV. RESONANCE BEHAVIOR 

Finally, we explore the possibility of the stochastic res- 
onance phenomena in the FFJJA under uniform oscil- 
lating currents. In numerical simulations, we compute 
the power spectrum P{u>) of the staggered magnetization 
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FIG. 5: Scaling plot of the dynamic order parameter versus 
the temperature for size L = 16(D), 24(0), an d 32(A) with 
Q,/2n = 0.08 and (a) I = 0.3; (b) I = 2.2. The data fit- 
ting has been made with the equilibrium values of the critical 
exponents (a) v — 0.82 and /3/u = 0.11; (b) v — 0.95 and 
P/v = 0.13. 



m{t) and its SNR, which is defined to be 
SNR= 101og 10 



(8) 



The signal S is given by the peak intensity of the power 
spectrum at the driving frequency and N represents 
the the background noise level, which is estimated by the 
average power spectrum around the signal peak. We can 
also define similar quantities for the higher harmonics of 
the driving frequency f2, the behaviors of which are found 
to be quite peculiar as described below. 

In Fig. ^, we display typical behavior of the power 
spectrum P(co) at low temperatures for several driving 
amplitudes. Note the substantial difference according to 
the driving current amplitude Iq, which may be under- 
stood in view of zero-temperature states. For Iq = 0.8, 



where there is no oscillation of the staggered magnetiza- 
tion m(t) at zero temperature, observed at low temper- 
atures are broad peaks around odd harmonics as well as 
sharp peaks at even harmonics, as shown in Fig. |](a). 
On the other hand, for higher values of Iq, correspond- 
ing to Q — at zero temperature, addition of thermal 
noise gradually suppresses the sharp peaks at odd har- 
monics while inducing broad peaks at even harmonics [see 
Fig. ^(b)]. When the system enters another ordered re- 
gion of < Q < 1 by further increase of To , the staggered 
magnetization exhibits only oscillations of even harmon- 
ics at zero temperature, which remain as sharp peaks at 
finite temperatures. Here, noise again helps to develop 
broad peaks around odd harmonics and interestingly, the 
third harmonics are dominant over the first harmonics 
[see Fig. ||(c)]. 

Such interesting behavior depending on the driving 
current amplitude is expected to bring about rich physics 
in the SR phenomena. The SNR computed at the first 
and second harmonics for 7 = 0.8 is plotted in Fig. |?] 
(a) and (b), respectively. It is demonstrated that both 
the broad peak at the first harmonics and the sharp one 
at the second harmonics, which shows up in the presence 
of the noise, exhibit the SR behavior. Remarkably, the 
resonance temperature for the first harmonics is found 
to be different from that for the second harmonics: The 
former is higher and the dynamic transition tempera- 
ture is located between the two resonance temperatures. 
For Io — 1.2, the SR behavior of the second harmon- 
ics is again observed as shown in Fig. ||(b). However, 
the SNR of the first harmonics decreases monotonically 
with the temperature [see Fig. ||(a)], indicating that the 
oscillations of the odd harmonics present in the zero- 
temperature state are simply suppressed by the noise. 
Similar argument explains the SR behaviors for Iq = 2.0, 
where only the even harmonics exist in the oscillations of 
the staggered magnetization at zero temperature. The 
thermal noise tends to reduce the oscillations of the sec- 
ond harmonics, which is manifested by the monotonic 
decrease of the SNR in Fig. ^(a). In contrast, the SNR 
of the third harmonics is initially enhanced by the weak 
thermal noise as shown in Fig. |](b). These interesting 
behaviors may be understood as follows: The thermal 
noise, which is not strong enough to destroy the chiral- 
ity lattice, is expected to help the chirality lattice move 
over one more lattice barrier if the threshold is not so 
large. Accordingly, for To = 1.2 corresponding to Q = 
at zero temperature, the even harmonics which are ab- 
sent in the ground state are induced by the thermal noise 
while the odd harmonics are generally suppressed. Via 
the same argument, one can expect that for Jo = 2.0 in 
the second ordered region, the odd harmonics exhibit the 
SR behavior with the even harmonics decreasing mono- 
tonically. At higher values of Io, similar SR behavior is 
expected for odd and even harmonics according to the 
zero-temperature oscillations of the staggered magneti- 
zation. 
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FIG. 6: Power spectrum of the staggered magnetization for 
size L = 16 and driving frequency Q/2tt = 0.08. The driving 
amplitude and the temperature are (a) I {) = 0.8 and T = 0.28; 
(b) 7 = 1.2 and T = 0.2; (c) I = 2.0 and T = 0.3. 
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FIG. 7: Signal-to-noise ratio versus the temperature in the 
system of size L = 16 with Q,/2n = 0.08 and 7 = 0.8. The 
signal-to-noise ratio is calculated (a) at the first harmonics 
and (b) at the second harmonics. Error bars have been esti- 
mated from standard deviations; lines are guides to eyes. 



V. CONCLUSION 

We have examined the dynamic properties of the 2D 
FFJJA driven uniformly by alternating currents, with at- 
tention paid to the dynamic transition and the SR phe- 
nomena. At low temperatures the chirality in the FFJJA 
displays antiferromagnetic ordering. To describe the dy- 
namic transition associated with the antiferromagnetic 
ordering of the chirality, we have introduced the dynamic 
order parameter, which is the staggered magnetization 
averaged over one period of the oscillating current, and 
investigated its behavior with respect to the amplitude 
and the frequency of the driving current as well as to 
the temperature. At zero temperature, as the driving 
amplitude is raised, the dynamically ordered phase and 
the disordered one have been observed to appear alter- 
natively. The saturation value of the dynamic order pa- 
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FIG. 8: Signal-to-noise ratio versus the temperature in the 
system of size L = 16 with Q/2tt = 0.08 and I = 1.2 The 
signal-to-noise ratio is calculated (a) at the first harmonics 
and (b) at the second harmonics. The leftmost data corre- 
spond to the temperature T = 0.001. 



FIG. 9: Signal-to-noise ratio versus the temperature in the 
system of size L = 16 with Q/2-k = 0.08 and I = 2.0 The 
signal-to-noise ratio is calculated (a) at the second harmonics 
and (b) at the third harmonics. 



rameter is unity in the first ordered region while it is less 
than unity in other ordered regions. It has also been ob- 
served that the increase of the driving frequency tends 
to shift the dynamically ordered regions toward larger 
values of the driving amplitude. Except for the scale, 
the overall behavior at zero temperature does not change 
qualitatively with the driving frequency. As the temper- 
ature is increased, the ordered phase undergoes a dy- 
namic phase transition to the disordered phase. We have 
obtained the phase diagram of the dynamic transition 
on the plane of the temperature and the driving ampli- 
tude. Through the use of the finite-size scaling analy- 
sis, the dynamic transition has been shown to belong to 
the same universality class as the equilibrium Z2 transi- 
tion in the FFXY model. We have also examined the 
SR phenomena in the system. To investigate the phe- 
nomena, we have calculated the power spectrum of the 



staggered magnetization, which shows characteristic be- 
havior according to the state at zero temperature. In 
the regions corresponding to finite values of the dynamic 
order parameter at zero temperature, broad peaks have 
been observed around odd harmonics as well as sharp 
peaks at even harmonics. On the other hand, in the re- 
gions with vanishing dynamic order parameters at zero 
temperature, the opposite situation has been observed. 
As manifested by the SNR at several harmonics, the odd 
and the even harmonics have been revealed to exhibit dif- 
ferent SR behaviors, which may be understood in view 
of the variation of the zero-temperature states with the 
driving amplitude. 
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